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ABSTRACT 

The discovery of quasi-periodic oscillations (QPOs) in low-mass X-ray bina- 
■ ries (LMXBs) has been reported and discussed in recent studies in theoretical and 

observational astrophysics. The Rossi X-Ray Timing Explorer has observed os- 
cillations in the X-ray flux of about 20 accreting neutron stars. These oscillations 

m 

to ~ 1200 Hz and correspond to dynamical time scales at radii of a few tens of 



O 



are very strong and remarkably coherent. Their frequencies range from ~ 10 Hz 



kilometers, and are possibly closely related to the Keplerian orbital frequency of 



matter at the inner disk. Almost all sources have also shown twin peaks QPOs in 
the kHz part of the X-ray spectrum, and both peaks change frequencies together. 

Bearing in mind these facts, several models explain the observed frequencies 
by the oscillation modes of an accretion disk. In the beat-frequency interpreta- 
+3 • tion, the upper kHz QPOs are associated with the Keplerian frequency of the 

accreting gas flowing around a weakly magnetized neutron star (B ~ 10 7 — 10 10 
G) at some preferred orbital radius around the star, while the lower kHz QPOs 
result from the beat between the upper kHz QPO and the spin frequency of the 
neutron star. In the relativistic precession model, while the upper kHz QPO is 
still the orbital frequency of matter at some radius in the disk, the lower kHz 
QPOs are considered as general relativistic precession modes of a free particle 
orbit at that radius. 

In this study we propose a new generic model for QPOs based on oscillation 
modes of neutron star magnetospheres. We argue that the interaction of the 
accretion disk with the magnetosphere can excite resonant shear Alfven waves 
in a region of enhanced density gradients. We demonstrate that depending on 
the distance of this enhanced density region from the star and the magnetic 
field strength, the frequency of the field line resonance can range from several 
Hz (weaker field, farther from star), to approximately kHz frequencies (stronger 
field, ~ 6 — 10 star radii from the star). We show that such oscillations are able 
to significantly modulate inflow of matter from the high density region toward 
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the star surface, and possibly produce the observed X-ray spectrum. In addition, 
we show that the observed 2 : 3 frequency ratio of QPOs is a natural result of 
our model. 



1. introduction 

Recent years have seen a surge in interest in studies of accreting neutron stars and 
black holes. One reason is that these objects provide a very unique window to a deeper 
understanding of the physics of strong gravity and dense matter. The complicated structure 
of these stellar objects, including strong gravity, gravit at ionally- captured plasma accretion, 
rapid rotation, bursting in accreted materials, thermal X-ray emission, and several other 
features makes a most fascinating object for both theoretical and observational investigations. 
One of the most interesting features of these objects is the interaction of accreting plasma 
with the stellar magnetic field. Matter is transferred from a normal donor star to the compact 
object. MHD models of the infalling plasma and the neutron star magnetosphere allow a 
relatively simple first approach to an understanding of plasma flow and the configuration of 
the magnetosphere around a neutron star. 

Research on plasma accretion by magnetic neutron stars began in the early 1970s with 
the discovery of periodic oscillations in X-ray fluxes of the stars Cen X-3 (Giacconi et al. 
1971; Schreier et al. 1972) and Her X-l (Tananbaum et al. 1972). Interest in these systems 
rose about two decades ago, after the first observation of quasi-periodic oscillations (QPOs) 
in the persistent X-ray fluxes of the luminous X-ray sources (van der Klis et al. 1985; Hasinger 
et al. 1986; Middleditch & Priedhorsky 1986). These oscillations are detected both in low 
and high frequencies with frequencies ranging from ~ 1 to ~ 100 Hz and from ~ 300 to 
~ 1200 Hz, respectively. The Rossi X-Ray Timing Explorer (RXTE) has observed persistent 
X-ray emissions about 20 accreting neutron stars in low- mass X-ray binaries (LMXBs). 
Generally, three different millisecond phenomena have now been observed in X-ray binaries: 
the twin kilohertz quasi-periodic oscillations (kHz QPOs), the burst oscillations, and the 
spin frequency of the accreting, low-magnetic field, neutron star. These kilohertz frequencies 
are the highest frequency oscillations ever seen in any astrophysical object, and are now 
widely interpreted as due to orbital motion in the inner accretion flow. The oscillations 
sometimes have one spectral peak. In most sources, however, the X-ray spectrum shows 
two simultaneous kHz peaks with a roughly constant peak separation, moving up and down 
together in time as a function of photon count rate. In this work, we focus on the Hz-kHz 
QPOs. 

The burst oscillations with frequencies in the 300 — 600 Hz range have been observed in 
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a number of the X-rays sources, usually during a type I X-ray burst due to a thermonuclear 
runaway in the accreted matter on the neutron star surface. Furthermore, the frequency 
of the burst oscillation increases by 1 — 2 Hz during the burst tail, converging in time to a 
stable frequency at the neutron star spin frequency. It is believed that the burst oscillations 
arise due to a hot spot or spots in an atmospheric layer of the neutron star. This has been 
explained by a decoupling of the atmospheric layer from the star and expansion by 5 — 50 m 
during the X-ray burst, while its angular momentum remains constant. Therefore, during 
the burst tail, the frequency increases due to spin-up of the atmosphere as it re-contracts. 
The asymptotic frequency corresponds to a fully re-contracted atmosphere and is the closest 
to the true neutron star spin frequency (see the recent review by Bildsten & Strohmayer 
(2003)). In a few sources, however, a rather large frequency drift, ~ 5 — 10 Hz, is observed. 
This drift is not consistent with the above model ( see Rezania & Karmand (2003) for more 
details). 

Although there are some detailed differences in different types of X-ray sources, the 
observed QPOs in all sources are remarkably similar, both in frequency and peak separation. 
(Six of the 20 known sources are originally identified as "Z" sources and the rest are known as 
"atoll" sources, see table 1. For more information about the atoll and Z sources, see van der 
Klis (2000).) Such a similarity shows that QPOs should depend on general characteristics 
of the X-ray sources which are common to all systems. In other words, the QPO can be 
regarded as a generic feature of the accreting neutron star. 

In the accretion disk of accreting binaries the material from the donor star moves around 
the compact star in near-Keplerian orbits. Depending on the binary separation, the disk 
radius varies from 10 5 to 10 7 km. Although there is still no certainty about the flow geometry 
in the inner emitting region, the extension of a part of the flow down into the emitting region 
in the form of a Keplerian disk is the basic assumption of most proposed models for accretion 
of matter onto low magnetic field neutron stars (see for example Miller et al. (1998)). Due to 
the interaction with the magnetic field of the star, radiation drag, and relativistic effects, the 
flow will be terminated at a radius rj n , the innermost stable orbital radius that is larger than 
radius of the compact star R s . Within r in matter may leave the disk and so the flow is no 
longer Keplerian. Remarkably, the characteristic dynamical time scale for materials moving 
near the compact object is comparable with the observed millisecond X-ray variability, ie. 
Tdyn = (r 3 /GM) 1/2 ~ 2 ms (r/lOOkm) 3 / 2 (1.4M /M) 1 / 2 . Such a natural periodicity is the 
foundation of most models for the observed QPOs. We will discuss these models in the next 
section. 

Another common feature in accreting binaries is the interaction of the accretion disk 
with neutron star's magnetosphere. In all sources, matters transfer from a normal donor star 
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to the compact object. This matter is accelerated by the gravitational pull of the compact 
object and hits the magnetosphere of the star with a sonic/supersonic speed. The MHD 
interaction of the infalling plasma with the neutron star magnetosphere, will alter not only 
the plasma flow toward the surface of the star, as assumed by current QPO models, but also 
the structure of the star's magnetosphere. The magnetic field of the neutron star is distorted 
inward by the infalling plasma of the Keplerian accretion flow. Since the gravity of the star 
confines the inward flow to a small solid angle ~ 1CT 2 ster, the magnetic field of the star 
will be more compressed in the disk plane than in other areas, see Fig. 1. Furthermore, 
in a more realistic picture, one would expect that the highly accelerated plasma due to the 
infall process would be able to penetrate through the magnetic field lines. Such materials 
will be trapped by the magnetic field lines and produce enhanced density regions within 
the magnetosphere. However, due to their initial velocities, they move along the field lines, 
finally hit the star surface at the magnetic poles and produce the observed X-ray fluxes (Fig. 
1). See Ghosh et al. (1977) for more detail. 

Besides modifying the geometry of magnetosphere, the compressional action of the ac- 
cretion flow can excite some perturbations in the enhanced density region. This can be 
understood by noting this fact that the inward motion of the accretion flow will be halted by 
the outward magnetic pressure at a certain distance from the star, the Alfven radius (Ghosh 
et al. 1977). In other words, the accretion flow pushes the stellar magnetic field toward 
the star until the inward pressure of the infalling plasma balances with the outward 

magnetic pressure Bp/8n. Here p and v r are the density and radial velocity of the infalling 
matter and B p is the poloidal magnetic field at the disk plane. Therefore, any instability in 
the disk at the Alfven radius would disturb such quasi-equilibrium configuration as well as 
the structure of magnetosphere. As an example, the interaction of the solar wind with earth 
magnetosphere excites resonant shear Alfven waves, or field line resonances (FLRs), along 
the magnetic field lines (Samson 1991). As a result, one might expect that such Alfven waves 
can be excited by the compressive accreting plasma in the magnetosphere of an accreting 
neutron star with the accretion flow playing the role of the solar wind. 

The plan of this paper is as follow: In section 2 we discuss some of the current models 
for the observed QPOs. Specifically, we will review the two models, the sonic-point beat 
frequency model and the relativistic precession model, in more detail. We outline an occur- 
rence of the resonant coupling in a magnetized plasma in section 3. Specifically, we review 
the excitation of shear Alfven waves by studying linear perturbations in MHD. Next, the 
occurrence of field line resonances (FLRs), as a result of a resonant coupling between the 
compressional and shear Alfven waves, is discussed. To illustrate the basic features of FLRs, 
the excitation of these resonances in a rectilinear magnetic field and in a dipolar field are 
considered. A dipolar topology is likely appropriate for the magnetosphere of a neutron star, 
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near the surface. Depending on the field line where the resonance occurs the eigenfrequency 
of the FLR is in range of several hundred Hz to kHz. In section 4, in order to consider a more 
realistic model for an accreting neutron star, we study an influence of the ambient flow along 
the field line on the excitation of FLRs. We demonstrate that in this case the eigenfrequency 
of the Alfven modes is modulated by the velocity of the plasma flow. Furthermore, in the 
presence of the field aligned flow, the plasma displacement parallel to the magnetic field lines 
is non-zero. These displacements, which are missing in case of zero ambient flow, might be 
responsible for modulating the motion of infalling materials toward the magnetic poles and 
producing the observed X-ray fluxes. A possible occurrence of more than one peak in the 
power spectrum is also discussed. In addition, we show that the observed 2 : 3 frequency 
ratio of QPOs is a natural result of our model. Section 5 is devoted to summarizing our 
results and further discussions. 



2. The QPO models 

Although the behavior of observed QPOs in the X-ray spectrum of accreting neutron 
stars in LMXBs are somewhat similar, a wide variety of models have been proposed for 
QPOs. Most, but not all of these models, involve orbital motion around the neutron star. 
These models include the beat-frequency models, the relativistic precession model, and the 
disk mode models. It is beyond the scope of the present work to discuss all these models. 
We point out, however, some of the main issues of the most prominent models. Further 
information can be found in a review by van der Klis (2000). 



2.1. The sonic-point beat-frequency model 

As discussed earlier, beat-frequency models are based on the orbital motion of the 
matter at some preferred radius in the disk. Alpar & Shaham (1985) firstly proposed a 
beat-frequency model to explain the low-frequency 1 horizontal branch oscillation (HBO) 
seen in Z sources (see also Lamb et al. (1985)). They used the Alven radius ta ~ 1-5 x 
10 cm M 17 2/1 ^(Mq/M)- 1 / 1 as the preferred radius 2 . Here M 17 is the mass of accretion 
rate in units of 10 17 g s _1 and /i 26 is the stellar magnetic field in units of 10 26 G cm 3 . 



1 Two different Low-frequency (< 100 Hz) QPOs were known in the Z sources, the 6 — 20 Hz so-called 
normal and flaring-branch oscillation (NBO; Middlcditch & Priedhorsky (1986)) and the f 5 — 60 Hz so-called 
horizontal branch oscillation (HBO; van der Klis et al. (1985)). 

2 Note that at ta we have pv% <~ Bp/4n. 
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Miller et al. (1996, 1998) proposed a beat-frequency model, the so-called the sonic-point 
beat-frequency model, based on a new preferred radius, the sonic radius. In this model, 
they assumed that near the neutron star there is a very narrow region of the disk in which 
the radial inflow velocity increases rapidly as radius decreases. Such a sharp transition in 
the radial velocity of plasma flow from subsonic to supersonic happens at the "sonic point" 
radius, r sonic 3 . This radius tends to be near to the innermost stable circular orbit (ISCO) 
4 j r isco; however, radiative stresses may change its location, as required by the observation 
that the kHz QPO frequencies vary. Comparing the HBO and kHz QPO frequencies, clearly 
^sonic <C r\, so part of the accreting matter must remain in near-Keplerian orbits well within 

In order to explain how such orbital frequencies can modulate the X-ray fluxes, they 
assumed that in the accreting gas some local density inhomogeneities "clumps" are created 
by magneto-turbulence, differential cooling of matters, and radiation drags. At r sonic these 
orbiting clumps gradually accrete onto the neutron star, following a fixed spiral-shaped 
trajectory in the frame corotating with their orbital motion and hit the star surface before 
they are completely sheared or dissipated by radiation drag. At the "footpoint" of a clump's 
spiral flow the matter hits the surface of the star and so the emission is enhanced. The 
clump's footpoint travels around the neutron star at the clump's orbital angular velocity, 
so a distant observer whose line of sight is inclined with respect to the disk axis sees a 
hot spot that is periodically occulted by the star with the Keplerian frequency at r sonic . 
This produces the upper kHz peak at u 2 . The high coherency of the QPO implies that all 
clumps are near one precise radius and live for several 0.01 to 0.1 s, and allows for relatively 
little fluctuations in the spiral flow. To explain twin peaks, it is also assumed that some 
of the accreting plasma is channelled by the magnetic field onto the magnetic poles that 
rotate with the star. The beat frequency at v\ occurs because a beam of X-rays generated 
by accretion onto the magnetic poles sweeps around at the neutron star spin frequency v s 
and hence irradiates the clumps at r sonic once per beat period, which modulates, at the beat 



3 Recently Zhang (2004) proposed a model for kHz QPOs based on MHD Alfven oscillations. He introduced 
a new preferred radius, 'quasi sonic point radius', where the Alfven velocity at this radius matches the 
Keplerian velocity. The author suggested that the upper and lower kHz frequencies are the MHD Alfven 
wave frequencies corresponding to different Alfven velocities (or different mass densities). Although his 
results have a good agreement with observations, no mechanism suggested in order to explain the excitation 
of Alfven wave oscillations through star's magnetosphere. Further, it is not clear how these Alfven wave 
frequencies modulate the X-ray flux coming from the surface of the star. 

4 In general relativity, no stable orbital motion is possible within the innermost stable circular orbit 
(ISCO), i?isco = GGM/c 2 w 12.5Mi.4m km. The frequency of orbital motion at the ISCO, the highest 
possible stable orbital frequency, is j^isco ~ (158O/Mi.4M ) Hz. 
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frequency, the rate at which the clumps provide matter to their spiral flows and consequently 
the emission from the footpoints. 

An important prediction of the sonic-point model is that Av = v 2 — V\ be constant at u s , 
which is contrary to observations (Lamb & Miller 2001). However, one obtains better consis- 
tency with observations by assuming that the orbits of clumps are gradually spiraling down. 
As a result, the observed beat frequency will be higher than the actual beat frequency at 
which beam and clumps interact. This can be understood by noting that during the clump's 
lifetime the inflow time of matter from clump to surface gradually diminishes. Therefore, 
the lower kHz peak gets closer to the upper one, and so Av decreases (Lamb & Miller 2001). 
Such a decrease is enhanced at higher X-ray luminosity L x because of a stronger radiation 
drag causes a faster spiraling-down, as observed (Lamb & Miller 2001). Due to complica- 
tions, it is hard to predict how exactly this affects the relation between frequencies which 
makes testing the model more difficult. 



2.2. The relativistic precession model 

Based on general relativity, it is well known that the motion of a free-particle in an 
inclined eccentric orbit around a spinning object experiences both relativistic periastron 
precession similar to Mercury's (Einstein 1915), and nodal precession (a wobble of the orbital 
plane) due to relativistic frame dragging (Thirring & Lense 1918). Noting these facts, Stella 
& Vietri (1998, 1999) proposed the relativistic precession model in which the high kHz QPO 
frequency v 2 is identified with the Keplerain frequency of an orbit in the disk (similar to the 
beat-frequency model) and the low kHz QPO frequency V\ with the periastron precession of 
that orbit. They also compared the nodal precession of the same orbit to the frequency of 
one of the observed low-frequency (v^p ~ 10 — 100 Hz) in LMXBs. 

In this model the kHz peak separation varies as Av = z/ 2 (l — 6GM /re 2 ) 1 / 2 caused by 
the relativistic nodal precession and the relativistic periastron precession (Stella & Vietri 
1998, 1999; Markovic & Lamb 1998). Here I is the star's moment of inertia and r the orbital 
radius. Morsink & Stella (1999) and Stella et al. (1999) noticed that stellar oblateness affects 
both precession rates and must be corrected for. Using acceptable neutron star parameters, 
they obtained an approximate match with the observed v 1 , v 2 and v LF relations provided 
v LF is twice (or perhaps sometimes four times) 5 the nodal precession frequency (Morsink & 
Stella 1999). However, to get a more precise match between model and observations, one 



5 It is suggested that such a relation could in principle arise from a warped disk geometry (Morsink & 
Stella 1999). 
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needs to use additional free parameters (Stella & Vietri 1999). 

In this model Av and v s are not expected to be equal as in a beat-frequency interpre- 
tations. An interesting result of the relativistic precession model is that Av should decrease 
not only when v 2 increases (as observed) but also when it sufficiently decreases. 

Although the relativistic precession model explains the three most prominent frequen- 
cies observed in the X-ray spectrum of LMXBs, with the three main general-relativistic 
frequencies characterizing a free-particle orbit, there are still some unanswered questions. 
For example, how stable are precessing and eccentric orbits in a disk, how is the flux modu- 
lated at the predicted frequencies, and why does the orbital frequency z/ 2 change with X-ray 
luminosity?. Furthermore, from this model relatively high neutron star masses (1.8 — 2M ), 
relatively stiff equations of state, and neutron star spin frequencies in the 300 — 900 Hz range 
are obtained. 



3. The magnetospheric model 

3.1. Field line resonances in the earth's magnetosphere 

In the last few decades there was a large expansion of the exploration of the near- 
Earth, plasma environment. One of the recent advances in the space plasma physics is the 
proof of the connection between FLRs and auroral dynamics (Samson et al. 1996). Samson 
ct al. (2003) use a combination of ground based and satellite observations with analytical 
and computational models to demonstrate a connection between FLRs and the formation 
of auroral arcs. Dobias et al. (2004) outline a possibility that the FLRs are a triggering 
mechanism for the nonlinear plasma instabilities that occur in the near-Earth magnetotail 
during magnetospheric substorms. Since FLRs are a generic feature exciting in magnetized 
plasmas with gradients in the Alfven velocity and with reflection boundaries, they are likely 
to occur in accreting neutron star magnetospheres as well. In the following sections we 
outline the occurrence of the resonant coupling within the ideal MHD approach and illustrate 
the resonance mechanism with some simple examples. We present results of computational 
models of shear Alfven waves in dipolar magnetic field, and we discuss the effect of an 
ambient flow parallel to the field. 



- 9- 



3.2. Magneto- hydrodynamic waves 

In this section we study MHD and applications to Alfven waves and FLRs. In general, 
the dynamics of a magnetized plasma is described by plasma density p = mn, plasma pressure 
p, gravitational potential $, velocity vector v and magnetic field B: 

dw/dt + v • Vv = -Vp/p + (l/47rp)B x (V x B) - V$, (la) 
dp/dt + V • (pv) = 0, (lb) 

An adiabatic equation of state with an adiabatic index 7, pj p 7 =const., is assumed in this 
paper in order to complete the above equations. 

In a steady state (ie. d/dt = 0), Eqs. (1) have been studied in detail regarding with the 
problem of stellar winds from rotating magnetic stars (Mestel 1961, 1968) and in connection 
with diffusing/flowing plasma into magnetospheres in accreting neutron stars (Eisner & Lamb 
1977, 1984; Ghosh et al. 1977). Obviously, in case of accreting neutron stars, v represents 
the inflow velocity of matter accreted to stars, while in stars with a stellar wind it represents 
the plasma outflow velocity. By decomposing the velocity and magnetic field vectors into 
poloidal and toroidal components: 

v = v p + B = Bp + B^, (2) 

one can obtain 

v p = (f/p)B p , (3) 

n = n a + (f/ P )(B <l> /w), (4) 

where / is the mass flux along a magnetic flux tube of unit flux and Q s = Q s ■ z is the 
angular velocity of the star (Ghosh et al. 1977). The subscripts p and <p denote poloidal and 
toroidal components, respectively, Q is the angular velocity of plasma at a distance w from 
the axis of the aligned rotator, and is a unit toroidal vector. We note that both / and Q 
are constant along a given field line. Equation (4) can be rewritten as 

n = n s + (v p /w)(B 4> /B p ), (5) 

where v p = v p ■ p is the magnitude of velocity along the poloidal magnetic field B p with 
magnitude of B p . Furthermore, Ghosh et al. (1977) have shown that at distances close to 
the Alfven radius, the poloidal component of the inflow velocity (for accreting systems) v p 
approaches the Alfven velocity va, ie. 6 

v 2 p (r A )=vi(r A ). (8) 



6 In order to estimate the poloidal component of the inflow velocity v p , one needs to integrate the mo- 
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The possible perturbations of a magnetized plasma and MHD waves are found by speci- 
fying the equilibrium configuration of the star and then solving Eqs. (1). This is a nontrivial 
problem, that is addressed by several investigators. However, in this paper we will be inter- 
ested in propagation of shear Alfven waves in the star's magnetosphere and their resulting 
FLRs. 



3.3. Field line resonances 

In MHD, waves can propagate in three different modes including shear Alfven waves, 
and the fast and slow compressional modes (Landau & Lifshitz 1992). In a homogeneous 
plasma, one can easily show that these three modes are linearly independent. However, in an 
inhomogeneous media these three modes can be coupled, yielding either a resonant coupling 
(Southwood 1974; Hasegawa 1976), or an unstable ballooning mode (Ohtani & Tamao 1993; 
Liu 1997). FLRs result from the coupling of the fast and the shear Alfven modes. 

The linear theory of the FLRs was developed by Chen & Hasegawa (1974) and South- 
wood (1974), and applied to auroral phenomena by Hasegawa (1976). Samson et al. (2003) 
developed a nonlinear model with a nonlocal conductivity model to explain the evolution 
of field aligned potential drops and electrical acceleration to form aurora. They studied the 
possible coupling between the fast compressional mode and the shear Alfven mode in an in- 
homogeneous plasma with radial gradient in the Alfven velocity va = B/^Anp. For further 
discussion of the field line resonances see for instance Stix (1992), and for the example of 
the numerical simulation see Rickard & Wright (1994). 

To illustrate the resonant coupling between the fast compressional mode and the shear 
Alfven mode we study the excitation of FLRs in an inhomogeneous plasma in two simple 
models for the magnetic field configuration, a rectilinear magnetic field model and a dipolar 



mentum Eq. (la): 

(l/2)(t>£ + n 2 vj 2 ) - n s tlm 2 - GM/r = const, (6) 

where M is the mass of the neutron star (Mestel 1968; Ghosh et al. 1977). In above equation the pressure 
term is neglected. Equation (6) shows conservation of energy in a corotating frame with the star, while in a 
nonrotating frame the extra term Sl s Qm 2 appears, that represents the work done by the magnetic field on the 
flowing plasma. However, as argued by Ghosh et al. (1977) the magnitude of v p inside of the magnetosphere 
is nearly equal to the free-fall velocity, ie. 



v p ~ (2GM/rfl 2 . 



(7) 
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magnetic field model. 



3.3.1. FLRs in a rectilinear magnetic field 

Plasma dynamics can often be described by the ideal MHD equations as: 

p(d/dt + v • V)v = -Vp + B x (V x B), (9a) 

dB/dt = V x (v x B), (9b) 

dp/dt + V • (pv) = 0, (9c) 

dipp-^/dt = 0. (9d) 

where we neglect the effect of gravitational attraction on the plasma here, see Eqs. (1) for 
more detail. Introducing Eulerian perturbations to the ambient quantities by 

p = Po + Sp, p = p + 5p, 
B = B + 5B, 5v = d£/dt, 

one can derive the linear wave equation by simplifying Eqs. (9) up to the first order in 
the perturbations. We choose coordinates such that the ambient magnetic field is in the z- 
direction. Without loss of generality, we assume that the gradients in the plasma parameters 
are in the x-direction. Setting the plasma displacement as 

£(r,i) = £(x)e~ l(ujt - k y y - k * z) , (10) 

Eqs. (9) reduce to (Harrold & Samson 1992) 

d 2 £ x F'(x)de x „, x> 

+ - F rr-r £ + G(a;)6 = 0, (11) 



dx 2 Fix) dx 



with 



and 



c 2 (c 2 - (v 2 + vi){k 2 y + A; 2 )) + A; 2 2 (A; 2 + fc» 2 

G{x) = Fi + viw - kiv\vi ' (12) 

2 _ U2 2 

where v\{x) = B 2 /4np and v$(x) = •yp/p, and ' = d/dx. Equation (11) yields two turning 
points at G(x) = (the compressional wave turns from propagating to evanescent), and two 
resonances at F(x) = 0. Close to the resonance positions, Eq. (11) can be approximated by 
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where xo is the position that the resonance occurs. 

In case of a strong magnetic field, or for a cold plasma (ie. Pfl U id/Pmagnetic C 1), one can 
put v s ~ 0. Then G{x) and F(x), (12) and (13), reduce to 



LO 



2 



G{x) = ^-K~K (15) 

and 



V A 



F{x) = : l( "^> 2 . (16) 
iv 2 - (fcjj + k 2 z )v{ 

In this case Eq. (11) has only one resonance (Alfven resonance) at 

o> 2 - k\v\ = 0, (17) 

which corresponds to the dispersion relation for the shear Alfven wave along the field line. 
At the resonance the incoming compressional wave reaches field line (or resonant magnetic 
shell) upon which the eigenfrequency of the standing shear Alfven wave (boundaries at the 
surface of the neutron star) is the same as the frequency of the incoming compressional, fast 
mode fluctuation. 

The FLR mechanism is generic and likely to occur in many astrophysical magneto- 
sphere. As a result, one would expect that the FLRs likely occure not only in the Earth's 
magnetosphere but also in the magnetospheres of accreting neutron stars. In the case of the 
Earth's magnetosphere, the source of energy for the resonant interaction is the interaction 
of the solar wind with the magnetosphere (Harrold & Samson 1992). In accreting neutron 
stars, the accreted plasma interacts with the stars' magnetosphere, allowing the compres- 
sional mode to propagate into the magnetosphere and flow along the field lines toward the 
magnetic poles. Such a compressional action of the accretion flow can excite resonant shear 
Alfven waves in the enhanced density regions filled by plasma flowing along the field lines. 
In section 4 we consider this mechanism more carefully to address its possible relation to 
those quasi periodic oscillations observed in accreting neutron stars in LMXBs. 



3.3.2. FLRs in a dipolar field 

In previous section we studied the excitation of FLRs in the simple box model (the rec- 
tilinear magnetic field). In this section, however, we consider a more realistic configuration, 
the dipolar magnetic field. Analytic solutions for FLRs in dipole fields can be found in Tay- 
lor & Walker (1984); Samson et al. (1996). Rickard & Wright (1994) outline computational 
(numerical) methods for modeling FLRs. 
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Let coordinates (xi, X2, x 3 ), such that the coordinate x± is directed along field lines, the 
coordinate X2 is in the direction of the radius of curvature, and the coordinate x 3 is in the 
azimuthal direction, see Fig. 2 7 . The equation of motion can be written in components as 



pdttCi 



0, 

B 
hih 2 

B 
hh 3 



d 2 (d 2 Bh 3 & + d 3 Bh 2 &) + d 1 



d 3 (d 2 Bh 3 Z 2 + d 3 Bh 2 S 3 ) + d 1 



h 2 



hih : 



-diBh 3 ^ 2 



fu 



-d L Bh^ 3 



(17) 
(17) 

(17) 



where hi, h 2 , h 3 are the metric (Lame) coefficients. We limited ourselves to a 2-dimensional 
problem in the X\X 2 (or xz) plane. Here d t = d/dt and <9, = d/dxi. 

In a strong magnetic field the usual definition of the Alfven velocity, v\ = B 2 /4np, is 
failed by resulting values larger than speed of light, c. To avoid of such confusion we define 
the relativistic Alfven velocity by using enthalpy 



h = 



c 2 p + B 2 /4tt 



(18) 



rather than the plasma density p. So the Alfven velocity shall be v\ = B 2 /4irh (Fig. 3). As 
shown in Fig. 3 close to the star the relativistic Alfven velocity saturates at the speed of 
light, ie. v\ ~ c. Such definition for the Alfven velocity is extensively used in studying of the 
relativistic jets dynamics in quasars. Figure 4 shows a computational model of the profile of 
the azimuthal component of the velocity. In a FLR with small azimuthal wavenumber, the 
displacement and velocity fields are predominantly azimuthal, with a 180° phase shift across 
the position of the resonance. 



4. Magnetohydrodynamics in the presence of ambient flow 

As discussed above, FLRs have been used to model auroral observations in the Earth's 
magnetosphere. Although the occurrence of these resonances is generic and likely to be 
excited in any magnetosphere with input of compressional of energy, one must carefully 
evaluate the differences between the Earth and neutron star magnetospheres. In the case of 
the Earth's magnetosphere, due to the small gravitational attraction of the Earth and also 
its large distance from the Sun, the solar wind more or less hits the whole Sunward side of the 



7 In dipolar coordinates, or coordinates close to dipolar, a notation (/i, v, <j>) is usually used where fi = 
cos#/r 2 , v = sin 2 0/r, and (r, 9, 4>) are spherical coordinates . 
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geomagnetosphere, and produces the so-called bow shock structure at the outer boundary. 
However, the strong gravity of the neutron star creates a supersonic converging flow long 
before the flow hits the star's magnetosphere. Such a localized flow is able to change the 
structure of the magnetosphere in local areas particularly in the equatorial plane. In addition, 
the highly variable nature of the exterior flow can change the magnetosphere's structure with 
time dramatically. The large flux of plasma stresses the outer star's magnetosphere and 
creates a relatively high plasma density throughout the magnetic field. The plasma then 
flows along the field lines, an interior flow, until it hits the star's surface near the magnetic 
poles, see Fig. 1. Furthermore, the strong magnetic field of the star, the rapid rotation, and 
the intense radiation pressure from the surface must be added to the problem. 

One of the most prominent differences between magnetospheres of accreting neutron 
stars and the Earth's magnetosphere is the presence of an ambient flow along the field 
lines where the resonance takes place (Ghosh & Lamb 1991). In this section we study the 
excitation of FLRs by considering such plasma flow in the magnetosphere. 

The presence of a flow v in the plasma adds more modes to the plasma waves. In 
general, such flow is a combination of plasma flow along the magnetic field lines, v p , and 
rotational motion of the plasma around the star with angular velocity f2, ie. v = v p + ft x r, 
see Eq. (2). The existence of v p itself allows a non-zero parallel component (relative to the 
magnetic field) of displacement that vanishes in a cold plasma with v p = 0. We note that this 
component is responsible for modulating the infalling plasma flow toward the star surface 
with FLRs' frequencies and producing the observed X-ray fluxes, see Eq. (22a) below. We 
will back to this point later. 

As we mentioned earlier, due to the relatively strong magnetic field in accreting neutron 
stars one needs to consider the relativistic Alfven velocity as 



See Eq. (18) for the definition of the enthalpy. Since the plasma pressure in neutron stars 
is negligible compared to magnetic pressure (Ghosh et al. 1977) we are using a cold plasma 
approximation in this section. Again we use a linear approximation to obtain a dispersion 
relation with resonant coupling. Since the gravity acts in a direction nearly perpendicular to 
field line for most of the infall of the matter, we can assume that the flow is nearly constant. 
This simplifies the derivation significantly. 

The linearized magnetohydrodynamic equations in the presence of an ambient flux can 
be obtained from Eqs. (9) as 




c 2 

c 2 p + B 2 /4n 47r ' 



(19) 



h 




-\- — — h v ■ V5v + 5v ■ Vv 
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+-!-5B x V x B + — B x V x 5B, (20a) 

47T 47T 

= V x (Sv x B + v x SB), (20b) 

at 
d 

( — + v • V)£p + 5v • Vp = -<y(5pV • v + pV ■ 5v), (20c) 

C6 

where <5v = d£/dt. The time scale of the ambient flow variations is much longer than time 
scale of the excited MHD perturbations, so dv/dt = 0. Furthermore, we consider the slow 
rotation approximation and so neglect the toroidal field B^, Eq. (2), to avoid complexi- 
ties. Such assumptions may not meet the actual configuration precisely. Nevertheless, these 
assumptions simplify our calculations significantly. 

Separating Eqs. (20) into parallel and perpendicular components relative to the ambient 
magnetic field and assuming perturbed quantities in the form of 



6(r,t) = ^r^e-^-^n) , (21) 



we find 



k iP' 

(-iu + iVn)£|| = T^ V • * - (£l * V >H > ( 22a ) 

uj — k\\v\\ 

(0- - *rf - *i)«x = - ^M^ft , (22b) 

(-iu + ik^SBw = iuB(V ± ■ £ ± ) + ■ V ± )B + (5B ± ■ V ± )v\\ , (22c) 

(-iu + ik\\v\\)SBj_ = uk\\B£ ± , (22d) 

(-iu + ik\\v\\)5p = iu(£ ± ■ V±)p + 7(-^||P£|| + iwpV± ■ £ ± ) , (22e) 

— r- = t—{< + v 2 A )V ■ £ p-£„ + €± • V ± V || , 22f 

n, w — || tu — K||f| (a; — ^||^||) 

where a; is the eigenfrequency, k\\ is the wave number in the parallel direction, v p = v\\B/B, 
v s — \fypjh is the sound velocity, and P = p + B 2 /8ir is the total plasma pressure (fluid + 
field) 8 . Note that we assumed that all ambient quantities are function of r_|_ only. Equation 
(22a) shows that the plasma displacement parallel to the ambient field does not vanish even 
in the cold plasma approximation (p = 0) 9 . The non-zero £y can affect and then modulate 
the motion of plasma along the field lines. Such modulation will occur at frequency of the 
resonant shear Alfven waves u as seen in the observed X-ray fluxes. 



8 In order to obtain Eq. (22f) we assumed that at equilibrium V(p + B 2 /8tt) « 0. This is valid for a very 
slowly rotating star. 

9 In case of zero ambient flow, v\\ =0, the parallel displacement £|| vanishes in the cold plasma 
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To analyze the problem analytically, we consider here again the rectilinear magnetic 
field configuration. We note that although this configuration may not be suitable for the 
accreting neutron star, it provides us a descriptive picture that can be applicable to QPOs. 
Setting the magnetic field in z-direction with gradients in the ambient parameters in the 
x-direction, f||(r_i_) = v p (x), and 

S(r,t) = 5{x)e- l{wt - k y y - k * z \ (23) 

Eqs. (22) in the cold plasma approximation (p = 0) can be combined into one differential 
equation in the form of Eq. (II) 10 

d% | 5'Qe)#» | (v - k z v p ) 2 - (k* + k*)vj 
dx 2 S(x) dx v\ 

where 

~( T ) - kzVpf - k 2 z vj 

" l ] u- k z v p (uj - k z v p ) 2 - (k 2 y + kl)v\ ■ { ) 

Note that here v p = v p (x) and va = va(x) are function of x. For v p — Eq. (25) reduces 
to one obtained in section 3.3, Eq. (13), as expected. The condition E(x) = yields the 
following modes: 

— k z (v p ± v A )- (26) 

The resulting modes are combination of a compressional mode k z v p due to the plasma flow 
and a shear Alfven mode k z VA- They produced two frequency peaks in the spectrum that 
vary with time if the velocities change with time. The value of k z ~ tt/L where L is the 
length of field line and is in order of radius of neutron star, L ~ R s . Therefore, for R s ~ 10 6 
cm and v a ~ v p ~ 3 x 10 8 cm s" 1 we get uj + ~ 1000 Hz that is comparable with frequencies 
of the observed quasi-periodic oscillations. In the limit v p va these modes become a single 
mode of a regular Alfven resonance. In the case of v a c these modes become consistent 
with a single mode (17). For a superalfvenic flow, ie. va <^ v p , however, the MHD Alfven 
waves and the resulting FLRs are suppressed due to propagation of the hydrodynamical 
wave, k z v p . Therefore, no FLRs are likely to occur at this level. This can be understood 
from Eq. (25) that for superalfvenic plasma flow reduces to 

v 2 k 2 

5(s)^ t *« (27) 
u - k z v p k 2 + kj 

where as u — > k z v p , E(x) — > oo. As a result, we expect that the FLR occurs where v p and va 
are comparable. Using the flow velocity and the shear Alfen velocity definitions (see below), 
one might expect that the FLRs likely occur at R s < r < 100i? s . 



10 



Note that in order to simplify our calculations we neglect the last term in the RHS of Eq. (22f). 
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Approximating the plasma inflow velocity with the free fall velocity v p ~ vg(r) and 
va ~ B(r)/y/A-Kp s where p ff = M / (v s Anr 2 ) is the free fall mass density, one can rewrite 
Eq. (25) as: 

w ± (r) ~ k z ( (2GM/r) 1/2 ± (B(r) 2 /47rp ff (r) ) 1/2 ) , 

~ k z (v p (R s )x-^ 2 ± ^(i?,)^- 9 / 4 ) , (28) 

where v p (R s ) = (2GM/R S ) = 1.6 x 10 10 (M/M Q ) 1 / 2 (i? s /10 km)^ 1 / 2 cm s" 1 and v A (R s ) = 
B/y/AnpsiRs) - 3 x lO^eM^^M/M©) 1 / 4 ^ 10 km)- 9 / 4 cm s" 1 are inflow and Alfven 
velocities at the surface of the star 11 . Here pL 26 is the magnetic field dipole moment at the 
surface of star in units of 10 27 G cm -3 , Myj is the mass of accretion rate in units of 10 17 g 
s -1 , x = r/R s , and R s is the radius of star. Figures 5 and 6 show the variation of purely 
Alfven resonance frequencies (v p = 0) or FLRs, as a function of distance from the star and 
magnitude of Alfven velocity. It is clear that the closer to the star and/or the larger the 
Alfven velocity the higher the frequency. Therefore, for an Alfven velocity say va ~ lc 
one can get a frequency about 1000 Hz at r ~ 2R S as we expected. Furthermore, since 
the Alfven velocity depends on the enthalpy (or density) of the plasma which varies from 
time to time due to several processes such as magneto-turbulence at boundaries during the 
accretion, the resulting FLRs frequencies will also vary with time. In addition, the position 
that the resonance takes place, is also subject to change in time due to the time varying 
accretion rate. Such time varying behavior causes time varying frequencies or quasi-periodic 
frequencies as observed. 

Setting oj~ as the observed peak separation frequency Au ~ u~, and u + as the upper 
QPO frequency u 2 ~ ou + , we can compare our results with observations. Figures 7 and 8 
show the variation of Au vs U\ and u 2 , respectively, along with the observed values. The 
curves are drawn by assuming the ratio Va/v p ~ .8 ar r = R s where from bottom to top 
correspond to upper frequency u 2 ~ 1170 Hz, 1350 Hz, and 1548 Hz at r = R s , respectively. 
As shown by observations, the value of Au decreases whenever the magnitude of u 2 decreases 
(sufficiently) and/or it increases. Such behavior is expected in our model (see Fig. 8). 



11 We note that to calculate the Alfven velocity at the surface of the star properly, one needs to use the 
relativistic Alfven velocity as defined in Eq. (19), ie. 

v A = i 17 =■ (29) 

yd + J u 2 Af- 1 (2GM) 1 /2 r -9/2/ c 2 

Using the above relation the Alfven velocity at the surface of the star will be ~ 0.8c, .997c, and c for 
/i = 10 26 , 10 27 , and 10 28 G cm 3 , respectively. In further distances from the star, however, the relativistic 
Alfven velocity reduces to its classical version. 
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Furthermore, the lower QPO frequency can be obtained from v\ — v 2 — Ais ~ 2k z v A . 
As a result, v\ to v 2 we obtain the frequency ratio 

v x jv 2 = 2v A /(v p + v A ) = 2/(1 + v p /v A ) . (30) 

It is clear that depending on the relation between the plasma flow velocity v p and the Alfven 
velocity v A , one can obtain different frequency ratios. For example, if the velocity inflow 
plasma is 2 (7/3 ~ 2.3) times bigger than the Alfven velocity, ie. v p = 2v A , one can get 
2 : 3 (3 : 5) ratio. For a free falling plasma along the lines of a dipolar magnetic field, 
v p (r) oc r _1 , whereas B 2 [r) / An p{r) oc r~ 9 / 2 , the condition v p ~ 2v A (v p ~ 2.3v A ) is satisfied 
at r ~ 3R S ~ 1.5r A . The occurrence of FLRs is very likely at this distance from the star. 

Such pairs of high frequencies with 2 : 3 frequencies ratio have been observed in the X-ray 
flux of almost all neutron stars in LMXBs 12 and black hole systems (the 3 : 5 ratio has seen 
in one source), and it would appear this feature is common to those systems. However, the 
existence of such rational ratios is still a mystery. In black hole systems, it has been suggested 
that such frequencies correspond to a trapped g-mode or c-mode of disk oscillation in the 
Kerr metric, see for example Kato (2001) and references therein. In neutron star systems they 
are explained as the fundamental and the first harmonic of the non-axisymmetric (m — 1) 
g-mode (Kato 2002, 2003). Further, Rezzolla et al. (2003a,b) studied small perturbations 
of an accretion torus orbiting close to the black hole and modelled the observed high QPO 
frequencies with basic p-modes of relativistic tori. They showed that these modes behave 
as sound waves trapped in the torus with eigenfrequencies appearing in sequence 2:3:4:... 
Abramowicz et al. (2003) also proposed that the observed rational ratios of frequencies may 
be due to the strong gravity of the compact object and a non-linear resonance between radial 
and vertical oscillations in accretion disks. 



5. Discussion 

In the present work, we have studied the interaction of an accretion disk with a neutron 
star magnetosphere in the LMXBs. The recent extensive observations reveal the existence 
of quasi-periodic oscillations in the X-ray fluxes of such stars. These oscillations, with 
frequencies range from 10 Hz to 1200 Hz, have been the subject of several theoretical and 
observational investigations. Based on theoretical models for the observed aurora in the 
earth magnetosphere, we have introduced a generic magnetospheric model for accretion 
disk-neutron star systems to address the occurrence and the behavior of the observed QPOs 



In Sco X-l the correlation line between two frequencies has a steeper slope than 2/3. 
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in those systems. In order to explain those QPOs consistently, we consider the interaction 
of the accreting plasma with the neutron star' magnetosphere. Due to the strong gravity 
of the star, a very steep and supersonic flow hits the magnetosphere boundary and deforms 
its structure drastically. Such a plasma flow can readily excite different MHD waves in the 
magnetosphere, such as the shear Alfven waves. 

In the Earth's magnetosphere, occurrence of aurora is a result of resonant coupling 
between the shear Alfven waves and compressional waves (produced by the solar wind). 
These resonances are known as FLRs that are reviewed in detail in section 3 of this paper. We 
argue that such resonant coupling is likely to occur in neutron star magnetospheres during 
its interaction with accreting plasma. We formulated an improved FLR by considering a 
plasma flow moving with velocity v p along the magnetic field lines. Such flow is likely 
to occur in neutron star magnetosphere (Ghosh et al. 1977). For a simple geometry, the 
rectilinear magnetic field, and in the presence of a plasma flow we found: (a) two resonant 
hydrodynamical-MHD modes with frequencies = k z (v p ±v a), (b) the resulting frequencies 
for k z ~ tt/R s ~ 3 x 10~ 6 cm -1 and typical flow and/or Alfven velocity ~ Ac will be in 
kHz range. Our results match the kHz oscillations observed in the X-ray fluxes in LMXBs. 
As shown in figures 5 and 6, the closer to the star and/or the larger the Alfven velocity the 
higher the frequency, (c) the quasi-periodicity of the observed oscillations can be understood 
by noting that due to several processes such as magneto-turbulence at boundaries and the 
time varying accretion rate, the FLR frequencies may vary with time, (d) a non-zero plasma 
displacement along the magnetic field lines £||. Such a displacement that oscillates with 
resulting frequencies u; ± modulates the flow of the plasma toward the surface of neutron 
star. As a result, the X-ray flux from the star will show these frequencies as well, (e) setting 
Av = uj~ and v 2 = uj + , one can well explain the behavior of the peak separation frequency 
Av relative to the upper QPO frequency v 2 . As observed, the value of Ais decreases as the 
magnitude of v 2 decreases and/or increases. Figure 8 clearly shows such behavior, (f) for 
v p ~ 2t>A at r ~ 3R S , the frequency ratio ui/u 2 is comparable with the observed frequency 
ratio 2 : 3. 

Interestingly, using the observed values of QPO frequencies, one can determine the 
average mass density and the magnetic field density of the star as 



k z v p (r) = (l/2)(i/ 2 + Av) = {y 2 - Vl /2), 
k z v A (r) = (l/2)(i/ 2 - Av) = v ± /2. 



(31a) 
(31b) 



Therefore, 



MJR] ~ {\/%ir 2 G)x\v 2 + Av) 2 , 

fi s ~ (327r A Gr l ^R 5 J 2 M 1 / 2 (M s /R 3 s )- 1 / A x l3 / A u 1 , 



(32a) 
(32b) 
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where k z ~ ir/r, /i s = B R^,x = r/R s , and B is the magnetic field strength at the surface 
of the star. In table 1, we calculate the average density M s /B? s and magnetic dipole moment 
Us of the star for a fixed value x = 10. Note that larger frequencies may occur at smaller 
distances, ie. x < 10. Our results are comparable with realistic neutron star parameters. 

Furthermore, our model is able to explain those low frequency (~ 10 Hz) quasi-periodic 
oscillations observed in the rapid burster such as MXB 1730-335 and GRO J1744-28 (Masetti 
et al. 2000). Such low frequencies can be extracted from the model by considering smaller 
inflow/ Alfven velocity and/or further distances from the star. 

Nevertheless, in order to avoid a number of complexities in our calculations, we used 
approximations such as slow rotation and cold plasma. These assumptions may put some 
restriction on the validity of our model and results. Future studies will be devoted to 
overcoming those restrictions. 

VR wishes to thank Mariano Mendez for kindly providing QPO data. This research was 
supported by the National Sciences and Engineering Research Council of Canada. 
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Table 1: Observed frequencies of kilohertz QPOs in Z and atoll sources 







v\ 


^2 


Al> 


M s /B? s 


Ms 


Source 


Type 


(Hz) 


(Hz) 


(Hz) 


(10 g cm A ) 


(10 G enr 3 ) 


ScoX-1 


Z 


565 


870 


307+5 


2.6 


3.8 


GX5-1 


z 


660 


890 


298+11 


2.7 


4.1 


GX 17+2 


z 


780 


1080 


294+8 


3.5 


4.8 


CygX-2 


z 


660 


1005 


346+29 


3.4 


4.1 


OA 34U+0 




565 


840 


339+8 


2.6 


3.8 


GX 349+2 


z 


710 


980 


266+13 


2.9 


4.5 


4U 0614+09 


atoll 


825 


1160 


312+2 


4.1 


4.8 


4U 1608-52 


atoll 


865 


1090 


225+12 


3.2 


1.7 


4U 1636-53 


atoll 


950 


1190 


251+4 


3.9 


5.4 


4U 1702-43 


atoll 


770 


1085 


315+11 


3.7 


4.8 


4U 1705-44 


atoll 


775 


1075 


298+11 


3.5 


4.8 


4U 1728-3 


atoll 


875 


1160 


349+2 


4.3 


5.1 


KS 1731-260 


atoll 


900 


1160 


260+10 


3.8 


5.4 


4U 1735-44 


atoll 


900 


1150 


249+15 


3.7 


5.4 


4U 1820-30 


atoll 


795 


1075 


278+11 


3.4 


4.8 


AqlX-1 


atoll 


930 


1040 


241+9 


3.1 


5.7 


4U 1915-05 


atoll 


655 


1005 


348+11 


3.4 


4.1 


XTE J2123-058 


atoll 


845 


1100 


255+14 


3.4 


5.1 



Low (ui) and high (1^2) QPO frequencies with corresponding peak separation (Az/) observed in Z and atoll sources. Only one 
series data is given as example, for further detail see van der Klis (2000). Values of M a / R? s and fi a arc calculated from Eq. (32) 
for a fixed value x = 10. We note that larger frequencies may be occurred at smaller distances, x < 10. 
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Fig. 1. — Side view of the magnetic field and accretion disk in accreting neutron stars. The 
neutron star's strong gravity causes a very high velocity flow toward the magnetosphere. 
As a result, the magnetosphere is pushed inward in the disk plane but balloons outward in 
direction away from the disk plane. Some of the plasma may leave the disk and flow along 
the field lines. 




Fig. 2. — Sketch of curvilinear coordinates for curvilinear or dipole magnetic field configu- 
rations. 




Fig. 3. — v\ in equatorial plane 
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Fig. 5. — Radial profile of the frequency of a FLR. The closer to the star the higher the 
frequency. Here i>a = -5c. 
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Fig. 6. — Frequency of a FLR at QR S from the star as function of Alfven wave velocity. The 
larger the Alfven velocity the higher the frequency. 
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Fig. 7. — Variation of peak separation frequency Au with the upper frequency v-i- The 
resulting curves are compared with the observed values for different system. The solid 
curves from bottom to top correspond to ~ 1170, 1350, and 1548 Hz at r = R s . The 
dashed curve corresponds to v-i = 1278 Hz at r = R s . 
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Fig. 8. — Variation of peak separation frequency Av with the upper frequency v^. The 
resulting curves are compared with the observed values for different system. The solid 
curves from bottom to top correspond to v-i ~ 1170, 1350, and 1548 Hz at r = R s . The 
dashed curve corresponds to v 2 = 1278 Hz at r = R s . 



